Extending the random-phase approximation for electronic correlation energies: The 
renormalized adiabatic local density approximation 
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The adiabatic connection fluctuation-dissipation theorem with the random phase approximation 
(RPA) has recently been applied with success to obtain correlation energies of a variety of chemical 
and solid state systems. The main merit of this approach is the improved description of dispersive 
forces while chemical bond strengths and absolute correlation energies are systematically underes- 
timated. In this work we extend the RPA by including a parameter-free renormalized version of 
the adiabatic local density (ALDA) exchange-correlation kernel. The renormalization consists of a 
(local) truncation of the ALDA kernel for wave vectors q > 2/cf, which is found to yield excellent 
results for the homogeneous electron gas. In addition, the kernel significantly improves both the 
absolute correlation energies and atomization energies of small molecules over RPA and ALDA. The 
renormalization can be straightforwardly applied to other adiabatic local kernels. 

PACS numbers: 31.15.E-, 31.15.ve, 31.15.vn, 71.15.Mb 
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Increasing computational resources has recently 
boosted a major interest in calculating electronic cor- 
relation energies from first principles using the adiabatic 
connection fluctuation-dissipation theorem (ACDF) [J- 
4]. The computational cost for such methods is much 
higher than traditional correlation functionals in density 
functional theory, but has the great advantage that it 
includes non-local effects and does not rely on error can- 
cellation between exchange and correlation. The RPA 
represents the simplest approach to ACDF calculations 
and has already been applied to broad range of elec- 
tronic structure problems (5l4ll|. While the non- locality 
of RPA makes it superior to semi-local functionals when 
dispersive interactions are important @, 0, EH 12-14|, 
the accuracy for molecular atomization energies is com- 
parable to that of the Perdew-Burke-Ernzerhof (PBE) 
functional @, @|. The relatively poor performance for 
atomization energies can be attributed to a deficient de- 
scription of short-range correlation effects. Furthermore, 
total correlation energies are severely underestimated in 
RPA and an accurate description of energy differences is 
highly dependent on detailed error cancellation. A sim- 
ple and intuitively appealing idea to remedy this problem 
was proposed by Yan et al. [TBI (RPA+), however, the 
method does not seem to improve atomization energies 
although total correlation energies are much better de- 
scribed 0, S 0| • From a perturbative point of view, RPA 
can be improved by including a screened second order 
exchange term (SOSEX) [lfj, which exactly cancels the 
one-electron self-correlation energy of RPA, albeit, with 
a significant increase in computational cost. In addition, 
it has been shown that RPA results can be improved by 
explicitly including single excitation terms, which correct 
the use of non-selfconsistent input orbitals 
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the description of short-range correlation effects by ex- 
tending RPA with an exchange-correlation kernel. For 
the homogeneous electron gas (HEG), this approach has 
been analyzed for a range of known exchange correlation 
kernels [if, EH and led to the construction of new adia- 
batic kernels fitted to reproduce the HEG correlation en- 
ergy [3, EH- S° far, it seems that for accurate total energy 
calculations, the non-locality of exchange-correlation ker- 
nels is very important, whereas the frequency dependence 
is less critical. Moreover, the pair-distribution function 
derived from any local approximation for the exchange- 
correlation kernel exhibits an unphysical divergence at 
the origin Q . While correlation energies are still well de- 
fined, the divergence makes it very hard to converge cor- 
relation energies based on local kernels. Recently, a fre- 
quency dependent exact exchange kernel has been shown 
to produce accurate correlation energies for atoms and 
molecules [2(|. However, the computational cost of this 
approach is significantly larger than that of RPA and the 
method may not be directly applicable to periodic sys- 
tems. 

In this letter we derive a non-local exchange- 
correlation kernel, which does not contain any fitted pa- 
rameters. The construction is based on a renormalization 
of the HEG correlation hole, which removes the diver- 
gence of the pair-distribution function and brings total 
correlation energies much closer to experimental values 
than both RPA and local approximations for the kernel. 

From the adiabatic connection and fluctuation- 
dissipation theorem, it follows that the correlation energy 
of an electronic system can be written 
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From the point of view of time-dependent density func- 
tional theory (TDDFT), it is natural to try to improve 



Here x KS is the exact Kohn-Sham response function and 
\ x is the interacting response function of a system where 
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the electron-electron interaction v has been replaced by 
Xv. Using TDDFT, one may express the interacting 
response function in terms of the Kohn-Sham response 
function as 



x x = x KS + x KS f&^<> 



(2) 



where fff xc = Xv + f£ c is the Hartree-exchange- 
correlation kernel at coupling strength A. The simplest 
approximation for fjj xc is the random phase approxima- 
tion where the exchange-correlation part is neglected. A 
natural next step is to include an adiabatic local approxi- 
mation for the exchange correlation kernel. In particular, 
one could try the ALDA kernel 



f ALDA [n]{v,v') = 5{r-v')f i 
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In the follow- 



where f ALDA [n] = £p(ne» c EG ) 

n— n(r) 

ing we will only consider the exchange part of the adi- 
abatic kernel, since it has the simplifying property that 
f£ = Xf x - Additionally, we expect the effect of includ- 
ing a kernel in Eq. (j2J will be dominated by the ex- 
change contributions. As it turns out, the kernel Eq. 
([3|) neither improves on total correlation energies [l| nor 
molecular atomization energies and is plagued by con- 
vergence problems related to the divergence of the pair- 
distribution function. 

For the homogeneous electron gas the problem is nat- 
urally analyzed in reciprocal space where an accurate 
parametrization of the correlation hole is known 0, [2l[ . 
In Fig. Q]we show the exact coupling constant averaged 
correlation hole of the homogeneous electron gas and 
compare with RPA and ALDAx results. Whereas RPA 
underestimates the value at a large range of g-values, 
ALDAx gives a reasonable description at small q but 
overestimates the value for q > 2k p. The ALDAx cor- 
relation hole becomes zero when /jy = 0, which hap- 
pens exactly at q = 2k p. The divergence of the pair- 
distribution function originates from the slowly decaying 
tail at large q where fgxc is complete dominated by the q- 
independent f ALDA 0. The full ALDA correlation hole 
is very similar to the ALDAx correlation hole displayed 
here % 

The correlation energy is essentially given by the 
integral of the coupling contant averaged correlation 
hole. Despite the divergent pair-distribution function, 
the ALDAjf correlation energy is well defined but con- 
verges very slowly due to the slow decay of g(q) at large 
q. From the shape of the correlation hole it is expected 
that RPA underestimates the correlation energy, while 
ALT) Ax overestimates it. Since the bad behavior of 
ALDAx primarily comes from large values of q it is now 
tempting to introduce a renormalized ALDAx correla- 
tion energy obtained by cutting the g-integral at the zero 
point of g(q). The result is shown in Fig. [5] along with 
RPA, PGG and ALDAx correlation energies. It is 
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FIG. 1: (color online). Fourier transform of the coupling con- 
stant averaged correlation hole for the homogeneous electrons 
gas. Left: r s — 1. Right: r s — 10. 
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FIG. 2: (color online). Correlation energy per electron of the 
homogeneous electron gas evaluated with different approxi- 
mations for f xc . 



seen that the renormalized ALDAx gives a remarkable 
improvement compared to RPA, ALDAjf and PGG. Ex- 
cept for the r s — > limit, it also performs better than 
the functionals proposed by Corradini et al 23] and the 
static version of the Richardson- Ashcroft kernel [24[ (not 
shown), which were fitted to quantum Monte Carlo data 
and derived from many-body perturbation theory respec- 
tively [H. 

For the homogeneous electron gas, the cutoff is equiv- 
alent to using the Hartree-exchange-corrclation kernel 



f r H Tc nA [n](q)=e[2k F -q)f^[n] 
Fourier transforming this expression yields 



/flic H ( r ) = f* A [n](r) +v r [n](r), 

fALDA 
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sin(2fci?r) — 2kpr cos(2kpr) 
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Since kp is related to the density, it is now straightfor- 
ward to generalize this to inhomogeneous systems. We 
simply take r — > |r — r'| and kp — > (37r 2 n(r, r')) 1 / 3 with 
n(r,r') = (n(r) + n(r'))/2. Thus, we obtain a non- 
local functional with no free parameters by performing 
a simple local renormalization of the correlation hole. It 
can be regarded as an ALDA functional where the delta 
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FIG. 3: (color online). Correlation energy of the valence elec- 
tron in Na evaluated with RPA, ALDA, and rALDA. The 
dashed lines show the values obtained with the functionals 
for the homogeneous electron gas using the average valence 
density of Na. 



function in Eq. §5§ has acquired a density dependent 
broadening. At large separation it reduces to the pure 
Coulomb interaction and it is expected to retain the ac- 
curate description of van der Waals interactions charac- 
teristic of RPA. For example, in a jellium with r s = 2.0 
two points separated by 5 A gives a renormalized inter- 
action v r [r s — 2](|r — r'|) = 0.97v(|r — r'|) and the mag- 
nitude of the Coulomb part of the kernel is a factor of 30 
larger than f^ ALDA . 

The renormalized ALDA functional has been imple- 
mented in the DFT code GPAW [HSl, which uses the 

The re- 



projector augmented wave (PAW) method (27 1. 
sponse function is calculated in a plane wave basis set as 
described in Ref . [28[ . The coupling constant integration 
is evaluated using 8 Gauss-Legendre points and the fre- 
quency integration is performed with 16 Gauss-Legendre 
points with the highest point situated at 800 eV. Since 
the kernel Eq. ([5]) is only invariant under simultaneous 
lattice translation in r and r', its plane wave representa- 
tion takes the form 



r rALDA 
JGG' 



(q) 



dr / dr'e- J ° r /(q;i-,r')e 



JG'-r 



(6) 



where G and G' are reciprocal lattice vectors, q belongs 
to the first Brillouin zone, and 



framework, it is not trivial to obtain the PAW corrections 
for a non-local functional and we use the bare ALDAx 
kernel to calculate contributions to the rALDA kernel 
from the augmentation spheres [28j j . 

As a first test of the functional for ab initio applica- 
tions, we have calculated the correlation energy of the 
valence electrons of bulk Na. We do not have a number 
for the exact value of the correlation energy but due to 
the delocalized nature of the valence electrons it is ex- 
pected that the result should be close to the correlation 
energy of the homogeneous electron gas at the average 
valence density of Na. This is supported by the close 
agreement between the RPA correlation energy of Na 
and the homogeneous electron gas [5^]. We found the 
rALDA calculations to be converged when two nearest 
unit cells were included. The result is shown in Fig. [3] 
as a function of plane wave cutoff energy along with the 
RPA and ALDAx results. As expected, RPA underesti- 
mates the correlation energy while ALDAx overestimates 
it. Again, one should note the slow convergence of the 
ALDAx calculation originating from the q-independent 
kernel. For plane wave implementations, an additional 
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problem is posed by the divergens of /, 
small densities. A particularly nice feature of the ker- 
nel (J5]) is that the small density divergence of ALDA is 
regulated. For example, for small r = |r — r' one obtains 



f'r 



ALDA 



(8) 



whereas ALDA diverges. 

The accuracy of molecular atomization energies by 
RPA is comparable to that of PBE, however, total cor- 
relation energies are typically severely underestimated. 
ALDA, on the other hand tend to overestimate total 
correlation energies by approximately the same amount. 
This is clearly seen for homogeneous systems displayed 
in Figs. [2] and [3] and the trend is also observed for inho- 
mogeneous systems. In Table H] we show a few examples 
of atomic and molecular correlation energies calculated 
with the rALDA functional and compared with LDA, 
PBE, RPA, and ALDA results. The ACDF correlation 
energies were calculated in a 6x6x7 A unit cell. The 
RPA and rALDA results were calculated at increasing 
cutoffs up to 400 eV and extrapolated to infinity. The 
ALDAx results were extrapolated from 1000 eV, but are 
still not well converged with respect to cutoff and repre- 
sent a lower bound on the absolute ALDAx correlation 



/(q;r,r') = l£ 



giq-Rij ig-iq-O- 



)/(r,r' + R ii ). (7) 



Here we have introduced the lattice point difference 



R, 



Ri — R, and the number of sampled unit cells 



N (k-points). /(q;r, r') is thus periodic in both r and 
r' and fGG' DA ( c should be converged by sampling a 
sufficient number of nearest neighbor unit cells. While 
the response function is calculated within the full PAW 





LDA 


PBE 


RPA 


ALDAx 


rALDA 


Exact 


H 


-14 


-4 


-13 


6 


-2 





H 2 


-59 


-27 


-51 


-16 


-28 


-26 


He 


-70 


-26 


-41 


-19 


-27 


-26 



TABLE I: Correlation energies of H, H2 and He evaluated 
with different functionals. Exact values are taken from Ref. 
[30t | - All number are in kcal/mol. 



4 





LDA 


PBE 


RPA@LDA 


RPAQPBE 


ALDA 


rALDA 


Exp. 


H 2 


-113 


-105 


-109 


-109 (109) 


-110 


-107 


-109 


N 2 


-268 


-244 


-225 


-224 (223) 


-229 


-226 


-228 


2 


-174 


-144 


-103 


-112 (113) 


-155 


-118 


-120 


CO 


-299 


-269 


-234 


-242 (244) 


-287 


-253 


-259 


F 2 


-78 


-53 


-13 


-30 (30) 


-74 


-39 


-38 


HF 


-161 


-142 


-122 


-130 (133) 


-157 


-136 


-141 


H 2 


-266 


-234 


-218 


-222 (223) 


-249 


-225 


-233 


MAE 


33 


10.1 


14.9 


8.4 


19 


3.7 





TABLE II: Atomization energies of diatomic molecules. The 
ALDA values are taken from Ref. 0] and experimental values 
(corrected for zero point vibrational energies) arc taken from 
Ref. [3l[ Results in brackets are from Ref. [(|. All number 
are in kcal/mol. The bottom line shows the mean absolute 
error for this small test set. 



energies. It is clear that the rALDA functional performs 
much better than both RPA and ALDA. 

The significantly improved total correlation energies 
are a very nice feature of the rALDA kernel. However, 
most physical properties depend on energy differences 
and the kernel is not of much use if it does not perform at 
least as well as RPA for such quantities. In Tab. HH we 
display the atomization energies of a few simple molecules 
calculated with different methods. The RPA@LDA and 
RPA@PBE columns show the sum of Hartree-Fock and 
RPA energies evaluated at self-consistent LDA and PBE 
orbitals respectively. Whereas the Hartree-Fock term is 
nearly independent of input orbitals the RPA correlation 
energies show a significant dependence on the ground 
state functional. This dependence is unfortunate since 
there is no obvious choice for the set of input orbitals. 
In contrast, when an adiabatic approximation for the 
exchange-correlation kernel is used, a consistent choice 
is the ground state functional from which the kernel was 
derived Q. In the present case of ALDA and rALDA 
we thus only consider calculations on top of the LDA 
ground state. For these molecules the rALDA kernel is 
seen to underbind by a few kcal/mol (F 2 excepted) but 
is superior to the RPA and ALDA results. 

The additional computational cost of calculating the 
kernel is insignificant compared to evaluating the non- 
interacting response function and inverting the Dyson 
equation. For a pure exchange kernel, it is possible to 
perform the coupling constant integration analytically, 
however, it involves an inversion of the non-interacting 
response function, which may become near singular at 
particular frequencies. The numerical coupling constant 
integration thus represents an additional computational 
cost compared to RPA calculations. 

In summary, we have presented a new parameter free 
exchange kernel for total correlation energy calculations 
within the ACDF formalism. The kernel largely cancels 



the self-correlation energy of RPA and seems to perform 
better than both RPA and ALDA for molecular atom- 
ization energies as well as for simple metals. Allthough 
more benchmarking is needed, these preliminary results 
indicate that the rALDA functional is clearly superior to 
RPA. In contrast to RPA, the functional has the very 
nice feature that it provides a consistent choice of input 
orbitals beyond the Hartree approximation. Finally, it 
will be straightforward to extend the kernel to include 
ALDA correlation, which might be expected to improve 
results further, but we will leave this to future work. In 
fact, the renormalization method naturally generalizes to 
all semi-local adiabatic approximations, which all suffer 
from the same pathological behavior in their pair distri- 
bution functions and the present work just represents a 
single example of an entire class of renormalizcd adia- 
batic exchange-correlation kernels. 

This work was supported by The Danish Council for 
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